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CHAPTER  1 
INTRODUCTION 

1.1  General  Introduction. 

The  aechanisa  of  the  foraatioa  of  a  detonation 
wave  in  a  coabustifcle  gas  mixture  has  been  studied 
extensively  in  the  past,  both  theoretically  and 
experinentally.  Whereas  aany  aspects  of  the  transition 
process  are  still  not  fully  understood,  the  process  of 
the  stable  detonation  wave  has  been  fully  clarified. 

The  transition  of  deflagration  to  detonation  is  a 
fairly  rapid  and  rather  complex  process.  Its  duration 
is  affected  in  part  by  the  teoperature,  pressure  ana 
composition  of  the  initial  reactants  and  by  both  the 
amount  of  turbulence  present  before  ignition  and  the 
amount  generated  by  the  flaae.  Another  important 
factor  is  the  shape  of  the  vessel  in  which  the 
combustion  tabes  place.  Burning  gas  aixtures  in  long 
cylindrical  tubes  is  a  coiaaon  method  for  investigating 
detonation  phenomena. 


From  profiles  of  the  average  wave  velocity  along 
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the  length  of  a  coabustion  tube,  the  detonation 
induction  distance  can  be  determined.  The  induction 
distance  is  ordinarily  defined  as  the  distance  froa  the 
ignitor  to  the  point  where  a  detonation  wave  forms. 

This  distance  varies  fron  one  gas  mixture  to  another, 
and  is  affected  by  the  initial  conditions. 

Hhen  a  coabustible  gas  mixture  contained  in  a 
cylindrical  tube  of  constant  diaaeter  is  ignited  at  one 
end,  a  flace  front  is  produced  which  propagates  into 
the  unburnt  aixture.  Because  of  the  temperature  rise 
behind  the  flame  front,  the  pressure  of  the  burnt  gas 
increases.  This  increase  in  pressure  leads  to  the 
formation  of  a  shock  wave  at  some  distance  ahead  of  the 
flame.  The  zone  between  the  front  of  the  shock  wave 
and  the  tail  of  the  combustion  wave  constitutes  the 
detonation  wave  during  the  transition  period. 

Behind  the  shock  wave,  the  unburned  gas  is 
compressed.  The  formation  of  the  shock  wave  by  the 
propagation  of  the  flame  front  is  similar  to  that 
produced  by  a  moving  piston. 

Because  of  the  heating  of  the  unburned  gas  behind 
the  shock  wave,  the  flame  speed  increases.  In  analogy 
to  the  piston,  the  expanding  coabustion  gas  does  work 
on  the  unburned  gas  behind  the  shock  wave  and,  thereby. 


produces  detonation  teaperatures  which  are  higher  than 
the  corresponding  adiabatic  flaae  teaperatures  at 
constant  voluae.  Eecause  of  this  teaperature  rise  the 
flaae  front  continually  accelerates  and  thus  produces 
additional  pressure  waves  which  increase  the  strength 
of  the  shock  wave.  At  some  point,  the  flaae  front  is 
travelling  fast  enough  to  catch  up  with  the  tail  of  the 
shock  wave.  The  resulting  structure  of  a  shock  wave 
and  adjacent  coabustion  wave  is  called  an  overdriven 
detonation  wave. 

In  this  state  the  detonation  wave  is  unstable 
because  it  has  not  attained  its  final  speed.  The  Hicfe 
nuaber  of  the  tail  of  this  detonation  wave  (with 
respect  to  the  burned  gas)  is  less  than  one.  The 
detonation  wave  is  considered  stable  after  the  tail 
Hach  number  reaches  unity.  The  wave  is  then  called  a 
Chapiran- Jouguet  detonation  vave. 

An  overshoot  in  coabustion  gas  pressure  and  shock 
wave  velocity  occurs  as  the  flaae  front  rushes  up  to 
the  tail  of  the  shock  wave.  The  aagnitudo  of  this 
pressure  overshoot  depends  on  the  voluae  of  the  heated 
but  unburned  gas  in  the  region  between  the  shock  wave 
and  the  flaae  front  at  the  tiae  of  aerging.  The 
aerging  process  is  fast  enough  that  the  unturned  gas  in 


this  region  explodes  almost  instantaneously. 

A  region  of  low  pressure  exists  at  the  ignitor  eni 
of  the  tube  as  a  result  of  the  forward  momentum 
imparted  to  the  combustion  gas  by  the  passage  of  the 
shock  wave.  This  low  pressure  region  propagates 
rarefaction  waves  upstreaa  which  catch  up  to  and  reduce 
the  strength  of  the  overdriven  detonation  wave.  These 
waves  reduce  the  speed  of  the  detonation  wave  by 
reducing  the  gas  pressures  within  it.  The  net  result 
is  a  continual  increase  in  the  Mach  number  of  the  tail 
of  the  detonation  wave. 

Eventually,  the  tail  Hach  number  becoaes  equal  to 
one,  and  the  expansion  waves  no  longer  weaken  the 
detonation  wave.  At  this  point,  the  detonation  wave  is 
considered  stable  and  proceeds  down  the  rest  of  the 
length  of  the  tube  at  a  constant  rate. 

1.2  Statement  Of  The  Problem. 

The  question  of  whether  induction  distance  is 
affected  more  by  changes  in  initial  pressure  or 
temperature  has  never  been  completely  resolved. 
Furthermore,  the  role  of  the  speed  of  sound  in  the 
unturned  gas  mixture  and  how  it  affects  the  length  of 


the  induction  distance  has  yet  to  be  fully  understood. 

A  study  is  undertaken  here  to  provide  better 
insight  into  the  relative  iaportance  of  the  initial 
pressure  and  temperature  on  the  induction  distance. 
This  is  done  oy  calculating  the  relative  euergy 
transfer  to  the  stable  detonation  wa7e  in  several 
hydrogen/oxygen/diluent  gas  aixtures.  The  induction 
distances  for  these  gas  mixtures  are  experimentally 
determined  by  burning  then  in  a  cylindrical  combustion 
tube  at  initial  conditions  of  room  temperature  and 
atmospheric  pressure. 

By  studying  the  initial  flame  front  Bach  number, 
more  insight  into  the  concept  of  detonation  induction 
due  to  flame  front  acceleration  is  obtained.  The 
initial  flame  front  Bach  number  is  calculated  for  the 
various  mixtures  in  this  study  based  on  experimentally 
determined  flame  speeds.  The  investigation  of  this 
parameter  results  in  a  conclusion  as  to  whether 
relatively  high  initial  Bach  numbers  correspond  to  low 
induction  distances. 
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There  is  such  literature  concerning  the  transition 
fro*  deflagration  to  detonation,  but  the  majority  of  it 
deals  with  the  process  occurring  in  solid  explosives 
rather  than  in  gaseous  aixtures  (references  1,  2  and  3 
for  exaaple) .  The  lain  literature  available  concerning 
transition  in  coabustible  gas  aixtures  still  have  not 
jet  resolved  the  guestion  of  what  aechaniss  causes  the 
detonation  induction  distance  to  change  (either 
lengthen  or  shorten)  as  the  initial  conditions  of  the 
aixture  are  changed. 

Experiments  have  verified  aany  times  that  the 
detonation  induction  distance  is  affected  by  changes  in 
initial  gas  aixture  pressure,  temperature,  composition, 
vessel  configuration  (i.e.  length  and  diameter)  and  the 
amount  of  turbulence  present  before  and  after  ignition. 
Early  work  by  Sokolik  and  Schelkin  (ref.  4)  showed  that 
initial  pressure  changes  the  length  of  the  induction 
distance.  These  investigators  performed  experiments 
using  hydrogen  and  oxygen  mixtures  as  well  as  various 


hydrocarbon  fuels  and  oxygen.  The  inltal  pressure 
ranged  from  30  to  300  aa  of  aercury  (.04  to  4 
atmospheres).  The  results  clearly  shoved  that 
induction  distance  decreases  as  initial  gas  mixture 
pressure  increases.  This  effect  of  initial  pressure 
has  been  confirmed  many  times  since  then  (references  8, 
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26,  27,  28,  29  and  30). 

Lafitte  (ref.  5)  did  various  early  experiments 
with  hydrogen/oxygen  mixtures  at  different  initial 
temperatures.  The  results  of  his  work  (which  are 
discussed  in  ref.  S)  showed  that  induction  distance 
increased  as  initial  temperature  increased.  The 
experiments  by  lafitte  involved  inital  temperatures 
ranging  from  15  to  350  d’egress  centigrade  (  288  to  623 
degrees  Kelvin). 

Edse  and  Lawrence  (ref.  6)  did  similar  experiments 
only  to  obtain  the  same  trend  in  induction  distance 
with  initial  temperature.  These  authors  used  initial 
temperatures  ranging  from  123  to  300  degrees  Kexvin. 


",  %**«%“  •„’> 

wC 


■  *  «  * .  * . 


V-V-V-V' 


•-.••v-v-' 


They  concluded  that  the  reason  tor  the  change  in  length 
of  the  induction  distance  in  a  combustible  gas  mixture 
depends  primarily  on  the  acceleration  of  the  flaxe 
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front  as  well  as  on  the  rate  at  which  the  pressure  of 
the  combusted  gases  increases.  This  conclusion  was 
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supported  by  use  of  an  empirical  equation  that 
inversely  related  the  detonation  induction  distance  to 
an  eapirical  factor  whose  dimensions  were  that  of  an 
acceleration.  It  was  therefore  postulated  that  this 
factor  was  sose  unknown  function  of  the  acceleration  of 
the  flaoe  front. 

Bollinger,  Fong,  Laughrey  and  Edse  (ref.  7) 
concluded  that  turbulence  is  an  important  aechanisa  by 
■easuring  the  detonation  induction  distance  in 
hydrogen/oxygen  aixtures  where  longitudinal  rod  inserts 
of  various  diaeeters  were  placed  in  the  coabustion  tube 
before  ignition.  Ihe  aixtures  were  ignited  at  a 
temperature  of  approxieately  313  degrees  Kelvin  and 
both  cod  diaeeter  and  initial  pressure  was  varied. 

They  further  concluded  that  as  both  initial  pressure 
and  cod  dianeter  increased,  detonation  induction 
decreased.  A  final  conclusion  was  that  turbulence 
increases  both  due  to  the  pressure  waves  generated  and 
the  resulting  interactions  as  the  coatustion  wave 
encounters  area  discontinuities. 

Bollinger  and  Edse  (ref.  8)  also  gave  another 
insight  into  the  effect  of  turbulence  on  the  detonation 
induction  distance.  This  conclusion  represents  a 
contradiction  to  soae  previous  views  which  is 


illustrative  of  the  problems  inherent  in  describing  the 
complicated  combustion  process  of  transition. 
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Briefly,  it  is  known  (from  the  work  by  Shchelkiu 
discussed  in  ref.  9)  that  as  tube  diameter  decreases, 
the  detonation  induction  distance  increases.  If  the 
amount  of  turbulence  per  unit  volume  of  gas  mixture  is 
considered,  this  ratio  decreases  as  tube  diameter 
increases.  Therefore,  it  should  be  expected  that 
detonation  induction  distance  should  increase.  This 
behaviour  would  be  expected  because  the  flame 
propagation  rate  decreases  as  the  amount  of  turbulence 
decreases  and  the  detonation  process  should  therefore 
occur  in  a  longer  period  of  time.  However,  as  stated 
above,  induction  distance  has  been  experimentally 
observed  to  decrease  as  tube  diameter  decreases. 

The  same  type  of  contradiction  is  obtained  when 
the  Beynolds  number  of  the  unburnt  gas  mixture  set  in 
motion  by  the  passage  of  the  shock  wave  is  considered. 
If  this  ratio  is  based  on  tube  diameter,  it  is 
proportional  to  the  diameter  and  decreases  as  tube 
diameter  decreases.  Since  the  Beynolds  number  is  in 
turn  proportional  to  the  amount  of  turbulence,  from 
this  reasoning  we  would  also  expect  the  induction 
distance  to  increase  as  tube  diameter  decreases.  Once 
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again,  the  expected  result  is  in  disagreement  with 
experimentally  observed  behaviour.  These  conclusions 
exemplify  the  need  for  a  better  understanding  of  the 
role  of  turbulence  in  the  induction  process. 

Lafitte  and  Shchelkin  have  shown  that  detonation 
wave  velocities  in  rough  tubes  are  only  40  to  50 
percent  that  of  those  observed  in  tubes  with  smooth 
walls.  Bybanin  (ref.  10)  suaaarizes  the  reason  for 
these  results  as  being  due  to  the  gas  mixture  being 
ignited  at  the  points  of  wave  reflection  from  the  wall 
roughness.  The  temperature  at  these  reflection  points 
has  been  shown  to  be  higher  than  the  temperature  behind 
the  shock  wave.  Combustion  propagates  toward  the 
center  of  the  tube  and  fills  the  cross  section.  Thus, 
the  combustion  occurs  before  the  mixture  has  time  to 
ignite  as  a  result  of  the  compression  process  behind 
the  shock  wave. 

The  question  of  whether  the  unbucct  but  shocked 
gases  between  the  shock  wave  and  the  flame  front  ignite 
spontaneously  or  by  a  less  rapid  process  whereby  the 
flame  front  merges  with  the  tail  of  the  shock  wave,  has 
been  investigated  quite  thoroughly.  Oppenheia,  Stern 
and  Ortiew  (ref.  11)  have  analyzed  in  detail  the 
results  of  experiments  by  Schmidt,  Steiaicktand  Neutert 


(ref.  12)  using  propane  and  air  mixtures  and  Schlexren 
photographs  of  the  coabustion  process.  The  idea  of 
pre-ignition  due  to  the  formation  of  high  teaperature 
areas  close  to  the  shock  wave  causing  ignition  of  the 
gas  aixture  between  the  shock  wave  and  the  flaae  front 
was  investigated.  Oppenheim,  Stern  and  Urtxew 
concluded  that  the  phenomenon  cf  pre-ignition  does 
cause,  but  is  not  neccessary  for  the  formation  of  a 
detonation  wave.  In  either  case,  it  was  further 
concluded,  an  overshoot  in  both  detonation  wave 
pressure  and  velocity  would  be  observed. 

Beyer,  Ortiev  and  Oppenheim  (ref.  13)  have 
concluded  that  the  gasdynamic  process  of  compression 
between  the  shock  wave  and  the  flame  front  contribute, 
at  most,  only  4  percent  to  the  enhancement  of  che 
transition  process  in  hydrogen/axr  mixtures.  Thus,  th 
phenomenon  of  transition  to  detonation  must  be 
associated  with  heat  or  mass  transfer  from  the  flaue 
front  although  the  effects  of  the  unsteady  boundary 
layer  should  also  be  consxdered.  These  authors' 
conclusions  are  based  on  the  results  of  the  calculated 
value  of  the  fractional  progress  of  the  induction 
process  which  takes  the  value  one  when  detonation 
occurs.  Values  that  they  calculated  for  this  fraction 
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were  at  nost  egual  to  .04  .  This  leans  that  a  particle 
in  the  shocked  but  unburnt  gas  nixture  would  have  only 
achieved  a  pressure  and  tesperature  of  4  percent  the 
value  that  it  has  at  detonation  if  the  transition 
process  were  the  result  of  gasdynaaic  coapressicn  only. 

Stehlow,  Crooker  and  Cusey  (ref.  14)  have 
concluded  by  both,  theory  and  ezperiaent  that  the 
accelerating  shock  wave  preceding  the  conbustion 
process  can  lead  to  the  foraation  of  *hot  spots'  (pre- 
ignition  points)  near  the  shock  wave  which  initiates  an 
explosion  and  produces  a  detonation  wave.  Both 
Strehlow  et  al.  and  Beyer  et  al.  conclude  that  the 
transition  process  can  be  nodelled  once  the  coaplicated 
kinetics  between  the  shock  wave  and  flaae  front  are 
fully  understood. 

Atkinson,  Bull  and  Shuff  (ref.  15)  have 
investigated  the  foraation  of  spherical  detonation 
waves  in  hydrogen/oxygen  aixtures.  They  introduced  a 
relation  for  calculating  the  detonation  induction 
distance,  [  ],  as  a  function  of  the  induction  tine,  t, 
and  the  velocity,  u,  of  particles  transaitting  tne 
conbustion  wave  just  prior  to  detonation,  [  ]  *  ut.  A 
plot  of  t  versus  f,  the  stoichioaetnc  ratio  of 
hydrogen  to  air  in  the  nixture,  was  nade  after  t  had 


been  calculted  foe  various  values  of  the  induction 
distance  obtained  from  measurements  and  u  was  assuaed 
to  equal  the  speed  of  the  stable  C-J  nave  for  each 
aizture.  Another  conclusion  was  that  the  aost 
detonable  aixture  of  hydroqen/air  occurred  when  #  uao 
approximately  equal  to  1.2  . 

Mori  in  the  area  concerning  the  development  of 
relations  for  the  the  speed  of  the  flaae  front  as  it 
accelerates  into  a  combustible  gas  aixture  has  not 
progressed  auch  further  than  analysis  of  the  flaae 
front  immediately  after  iguition.  Studies  by  dr tie*, 
Laderaan  and  oppenheia  (ref.  16)  have  produced 
relations  that  are  based  on  the  principal  that  the 
particle  velocity  of  the  shocked  but  unburnt  gases, 
prior  to  when  the  flame  front  reaches  the  *alls  ot  the 
detonation  tube,  is  a  linear  function  of  the  flaaefcont 
area.  In  their  investigation,  these  authors  have 
considered  cases  in  which  the  flame  front  initially 
propagates  away  from  the  ignitor  in  a  hemispherical 
shape  and  where  it  propagates  in  the  shape  of  a 
prolate  spheroid.  The  justification  for  using  such 
shapes  is  that  they  are  based  on  photographic  records 
of  the  coabustiou  process.  Use  of  their  method  is 
dependent  on  the  experimentally  determined  values  of 


the  initial  area  of  the  flaae  front  and  its  initial 
average  velocity,  lhat  is  iaportant  about  this 
research  is  the  emphasis  on  the  importance  of  the 
initial  Bach  number  of  the  flame  front. 

Other  work  by  Edse  and  Strauss  (references  17,  Itt, 
19  and  20)  has  produced  data  concerning  the  flame 
speeds  of  various  carbon  aouoxide/air  mixtures  at  hijh 
ambient  pressures.  The  data  obtained  by  these  authors 
has  shown  that  flaae  speeds  for  carbon  aonoxide/air 
mixtures  decrease  with  pressure,  whereas  the  flame 
speeds  of  hydrogen/oxygen  mixtures  increase  with 
pressure.  A  major  problem  in  experimentally 
determining  the  flaae  speeds  of  mixtures  that  behave 
like  the  carbon  aouoxide/air  system  is  controlling 
flashback.  In  other  words,  higher  ambient  pressures 
reguire  higher  gas  speeds  which,  coupled  with  the  lower 
flame  speeds,  can  make  it  impossible  to  maintain  a 
flaae.  Empirically  derived  relations  for  the  flame 
speed  can  therefore  not  be  applied  to  pressures  and 
temperatures  that  are  much  higher  than  those  at  which 
the  experiments  were  conducted. 

Since  it  is  possible  to  apply  the  energy, 
continuity  and  Hugoniot  equations  to  the  combustion 


process  within  a  stable  detonation  wave,  sever ji 


■ethods  exist  in  the  literature  for  calculation  of  C-J 
detonation  wave  paraaeters.  An  early  aetbod  used  by 
Levis  and  Friauf  (ref.  21)  calculates  values  based  on 
the  frozen  speed  of  sound  in  the  coabusted  gases  at  the 
tail  of  the  detonation  vave.  This  aethod,  which  is 
generally  considered  to  be  the  first  developed,  does 
not  include  the  presence  of  atonic  oxygen  in  the 
coabusted  gases.  Ihe  authors  believed  that  the  amount 
of  atoaic  oxygen  forned  was  not  enough  to  affect  the 
results  of  their  calculations,  however,  if  one  wishes 
to  obtain  a  coaplete  description  of  the  state  of  the 
coabusted  gases,  atoaic  oxygen  should  be  taken  in  to 
account  since  in  aixtures  involving  only  hydrogeu  ana 
oxygen,  the  sole  fraction  of  atonic  oxygen  is  usually 
on  the  order  of  .02  which  is  not  negligible. 

Berets,  Greene  and  Kistiakovsky  (ref.  22)  and  Dunn 
and  Holfsou  (ref.  23)  base  their  methods  on  the  frozen 
speed  of  sound  in  the  coabusted  gases  and  therefore 
obtain  values  of  the  teapecature  at  the  tail  of  the 
detonation  wave  which  is  less  than  actual.  Coaparison 
of  teaperatures  calculated  by  the  above  authors  has 
shown  differences  as  auch  as  130  degrees  Kelvin.  Ihe 
reason  for  these  differences  is  generally  attributed  to 
their  use  of  theraodynaaic  functions  which  ace  not  of 


the  saae  accuracy. 

Eisen,  Gross  and  Bivlin  (ref.  24)  iotroduced  a 
■ethod  that  bases  results  on  toth  the  frozen  and 
equilibrium  speed  of  sound  in  the  coabusted  gases. 
Bollinger  and  Edse  (ref.  25)  introduced  a  simplified 
■ethod  for  calculating  stable  detonation  wave 
parameters  in  hydrogen/oxygen  mixtures  in  which  tu.s 
Hugoniot  equation  for  a  cheaically  reacting  gas  aixture 
is  solved  for  the  ratio  of  the  pressure  of  the  uuburnt 
to  burnt  gases  at  the  tail  of  the  detonation  wave. 

This  iterative  procedure  calculates  the  partial 
pressures  of  each  specie  present  in  the  coabustion  gas 
and  produces  correct  results  which  are  obtained  without 
having  to  calculate  either  the  frozen  or  equilibrium 
speed  of  sound. 

The  procedure  presented  in  chapter  three  of  this 
paper  iterates  until  the  C-J  condition  at  the  tail  of 
the  detonation  wave  is  obtained.  Equations  for 
calculation  of  paraaeters  in  gases  containing  up  to 
eleven  species  is  presented. 


CHAPiEH  THREE 


THEORETICAL  ANALYSIS 


This  chapter  deals  with  the  calculation  of  stable 
detonation  wave  parameters  in  combustible  gas  mixtures 
An  iterative  method  is  presented  for  determining  the 
pressure  and  temperature  of  a  combustion  gas  at  the 
tail  of  a  stable  detonation  wave,  and  the  associated 
wave  velocities  and  gas  speeds.  Eguations  for 
calculating  the  combustion  gas  composition  in  several 
hydrogen/oxygen/diluent  mixtures  are  presented  where 
the  diluent  gas  is  either  nitrogen,  carton  dioxide  or 
some  iuert  species. 

3.1  General  Method. 

The  Hugoniot  eguations,  the  energy  and  continuity 
eguations,  and  the  equation  of  state  are  useful  in 
describing  the  kinetic  and  gasdynamic  processes  that 
occur  as  a  stable  detonation  wave  passes  through  a 
combustible  gas  mixture.  The  fact  that  we  calculate 


parameters  for  statle  detonation  waves  enables  us  to 
check  the  accuracy  of  tne  results  since  we  know  that 
the  velocity  of  the  tail  of  a  stable  wave  is  equal  to 
the  equilibrium  speed  of  sound  in  the  gases  behind  it. 

Since  the  thermodynamic  functions  of.  the 
combustion  gas  are  complex  functions  of  temperature, 
and  the  composition  is  uependent  on  both  temperature 
and  pressure,  an  iterative  procedure  is  required.  Tne 
procedure  begins  by  specifying  initial  estimates  of  the 
temperature  and  pressure  of  the  combustion  gas.  These 
estimates  are  continually  improved  until  the  calculated 
value  of  the  tail  Bach  number  is  unity.  Also,  it  is 
assumed  that  the  combustion  process  occurs  in  both 
thermal  and  chemical  equilibrium. 

Since  this  study  deals  with  hydrogen/oxygen  based 
mixtures  that  have  a  third  gas  additive  of  either 
carbon  dioxide,  nitrogen  or  an  meet  gas,  the  following 


equilibria  are  considered: 
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In  the  calculation  of  the  overall  sole  nuwber 
change, 

Av(1)  -Zv^  (3.1-1) 

where  sole  numbers  on  the  left  side  of  the  reactions 
are  considered  negative,  and  the  sole  nunters  ou  the 
right  side  are  considered  positive.  Also,  i  suns  over 
all  of  the  species  in  the  particular  reaction. 

For  instance,  if  the  gas  aixture  in  consideration 
was  +  +  sight  order  the  equilibria  in 

the  following  wanner. 
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numbering  the  equilibria  will  wake  the*  easier  to 
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reference  while  calculating  the  elements  of  several 
subsequent  aatrices. 


in  initial  estiaate  of  the  pressure  at  the  tail  of 

the  detonation  wave  can  he  obtained  by  considering  the 

effect  of  heat  addition  to  a  supersonic  flow  of  a 

calorically  perfect  gas.  in  particular,  when  y  Joules 

of  heat  ace  added  to  one  kilogcaa  of  the  moving  gas, 

# 

the  pressure  ratio  across  the  resulting  theraal  wave 
is. 
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(3.1-2) 


The  dimensionless  heat  release  factor,  r" t~  » 

P  1 

can  be  expressed  in  terms  of  the  combustion  enthalpy  of 
the  fuel  and  the  dimensionless  specific  heat  of  the 
combustion  gas  as  follows. 
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where 


|^Hcomb,f  | 


is  the  absolute  value  of  the 


combustion  enthalpy  of  one  sole  of  the  fuel  at  the 


unbucnt 


gas  teapecature  Ti  .  The  v.  are  the  mole 

1  •  o 


numbers  of  the  constituents  of  the  uudissocxated 
combustion  gas  obtained  from  the  stoichiometric 

J  y  V  <■ 


equation  for  the  reaction,  yj[Ji  iS  the  dimensionless 
specific  heat  of  species  i  between  the  teaperatures 
and  .  The  empirical  correctxon  factor  6  takes 

dissociation  into  consideration  and  is  equal  to 
approximately  .5  .  A  reasonable  estimate  for  7  is 
approximately  1.3  . 

Next,  the  specific  volume  of  the  combustion  gas 
can  be  estimated  from  the  following  relation. 
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Finally,  the  initial  estimate  of  the  comLusted  gas 
temperature  is  obtained. 
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Id  equations  (3- 1-2)  to  (3. 1-5),  the  subscript  1 
refers  to  the  conditions  of  the  unburnt,  unshocked  <jjs 
sixture.  The  subscript  3  refers  to  the  conditions  at 
the  tail  of  the  detonatiou  wave.  In  subsequent 
equations,  the  superscript  ESI  refers  to  an  estimated 
value  of  a  quantity  which  will  equal  the  actual  value 
at  the  end  of  the  iteration. 

EST 

With  these  initial  estiaates  of  temperature,  , 

EST 

and  pressure,  P-j  ,  the  equilibriums  constants  tor  the 

dissociation  equilibria  are  calculated.  Table  (3-1) 

lists  the  K-relations  for  the  equilibria  of  page  18  . 

T 

The  is  the  coefficient  of  the  equilibrium  constant 
of  species  i  at  temperature  T  and  is  interpolated  iron 
Tables  (3-2)  and  (3-3). 

Following  calculation  of  the  equilibriums 
constants,  the  aole  fractions  of  the  species  in  the 
coatustion  gas  are  calculated.  Since  the  equations  tor 
calculating  the  composition  ace  unique  to  the  auount 
and  types  of  species  present  in  the  coabustion  gas  tor 
particular  initial  gas  mixtures,  it  will  suffice  at  the 
present  time  to  refer  to  the  composition  equations  of 
section  3.2  for  the  aixtures  considered  in  this  stuoy. 
These  coaposition  calculation  procedures  require  at 
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least  an  initial  estiaate  of  the  aole  fraction  of 
diatonic  oxygen  in  the  coabustion  gas.  Gas  aixtures 
containing  nitrogen  and  carbon  (as  well  as  hydrogen  and 
oxygen)  reguire  initial  estimates  of  both  diatonic 
oxygen  and  diatonic  nitrogen. 

Once  the  coaposition  for  the  estinated  teaperature 
and  pressure  has  been  deternined,  the  molecular  mass  of 
the  coabustion  gas  is  calculated, 

1ST 

97? 3  -  ZnL97?i  (3.1-6) 
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Followed  by  the  reduced  enthalpy. 


\ 

test 

* 

'CG 

o 

l3  *©1^(  <3,‘ 


where  the  individual 


(»). 


are  interpolated  fcra 


Tables  (3-4)  to  (3-5).  The  universal  gas  constant,  <6,  , 

is  egual  to  8314.33  J/Kaol  K,  and  the  specific  gas 

constant  is  egual  to, 
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The  absolute  foraation  enthalpies 


# 


are 


obtained  froa  Table  (3*6). 


Froa  the  equation  of  state,  we  deteroine  the  ratio 
of  the  specific  voluae  of  the  coabustion  gas  to  that  of 
the  unburnt,  unshocktd  gas  aixture  to  be, 

(*1  p3ST<^?CG 


Froa  the  Hugoniot  equation  for  a  chenically 
reacting  gas,  we  calculate  a  value  for  the  pressure  at 
the  tail  of  the  detonation  wave. 


X\  *  (1+"^)]  (3‘1'10) 


If  the  absolute  value  of  the  difference  between 
the  estiaated  and  calculated  pressures  is  large 
(greater  than  .0000001  for  instance),  we  ^eestiaate  the 
pressure  in  the  following  Banner, 
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where  x  is  approxiaately  eguai  to  .5. 


The  coapositioa  is  recalculated  using  the  net* 


estimate  of  pressure  and  the  old  estiaate  of 
teaperature.  If  the  difference  is  snail,  we  proceed  by 
calculating  the  dioensionless  specific  heat  of  the 
coabustion  gas. 


(3.1-12) 


are 


where  the  individual 
Tables  (3-7)  to  (3-8)  . 

The  frozen  ratio  of  specific 


interpolated  froa 


heats  is  calculated 


(3-1-13) 


Since  the  OJ  condition  is  based  on  the 
eguilibriua  speed  of  sound  in  the  coabustion  gas  at  the 
tail  of  the  detonation  wave,  the  effective  ratio  of 
specific  heats  aust  be  calculated.  This  ratio  is  a 
function  of  the  specific  heat  froa  eguation  (3.1-12),  a 
cheaically  reacting  specific  heat  ratio  and  the  effect 

of  the  local  sole  nuaber  changes  as  sound  waves 

\ 


and 


The  inverted  I 


square  matrices  |*i,j 


[bl.j]  *  are  defined  as  the  inversions  of  the 


•atrices  whose  eleaents  are. 


r  t  ^  vi1}  • 

M  -  E  1 


(3-1-15) 


[*1.J  *  [bl.j]  -  Avi1)#*vij)  (3.1-16) 


The  letters  1  and  j  refer  to  the  chemical  change. 
For  exaaple,  if  the  equilibria  on  page  21  are  being 
considered. 


1  -  OH,0tHtC02,H20  j  -  0H,0,H,C02,H20 

Since  the  reactions  are  nuabered,  the  subscripts  1  and 
j  can  be  written  as  numbers,  for  example. 


bl,l"bOH.OH  1  bl,2"bOH.O  1  bl,3"b0H,H  ,tc *  *  • 


The  elements  of  the  row  aatrix 
I  Ah13)i 


and  column 


aatrix  If  I  ,  are  obtained  froa  Table  (3-9) 


The  values  of  the  row  and  column  aatrices. 


AH<» 


-  Av' 


are  easily  obtained  since  the  “ jj^ — and  Avu;  ace  known 
Likewise,  the  coluan  natrix  |av^^J  is  easily 
deteroined. 

Once  the  effective  ratio  of  specific  heats  is 
calculated,  the  equilibciua  speed  of  sound  is  obtained 


and  the  Hugoniot  equation  for  a  nocaal  shock  wave 
produces  the  velocity  of  the  bead  of  the  detonation 
wave  as  follows. 


(3-1-18) 


The  continuity  equation  produces  the  velocity  of  the 
tail  of  the  detonation  wave. 


w 


3  * 


(3-1-19) 


If  the  absolute  value  of  the  difference  between  w-j 
and  wa,3  is  large  (greater  than  .0000001  for  instance) , 
we  Bake  a  new  estinate  of  the  teaperature  frost  the 
following  relation, 

fEST(n+l)  .  jtSKn)  .  (  »3  .  Vj)  (3.J-M) 

where  x™,  is  approximately  egual  to  .5  -  The  entire 

A3 

iteration  routine  is  begun  again  until  the  magnitude  of 
the  difference  between  and  ^  is  very  snail. 

Lastly,  the  velocity  imparted  to  the  combustion 
gas  by  the  passage  of  the  shock  wave  is  obtained  as 
follow* , 

u3  *  *i  "  wy  (3-1-21) 

The  following  section  3.2  contains  the  composition 
eauati-ns  that  are  used  to  compute  the  mole  fractions 
of  the  combustion  qas.  They  are  presented  in  an 
algorithm  form  to  demonstrate  the  iterative  nature  of 


the  calculation 


3.2  Coo  position  Equations 


To  introduce  the  equations  for  calculating  the 
coabustion  gas  coaposition,  the  case  of  the 
hydrogen/oxygen/carbon  dioxide  aixture  will  be 
considered.  First,  the  ratio  of  the  global  sole  number 
of  diatoaic  oxygen  to  that  of  diatoaic  hydrogen  is 
calculated. 


(3-2-1) 


This  ratio  remains  constant  during  the  coabustion 
process.  It  is  calculated  during  each  iteration  and  is 

cmlc 

denoted  by  \  .  One  of  the  conditions  for  ending  the 

composition  iteration  is  the  requirement  that, 

I  calc  I 

f°2  '  SI  <  ‘°000001 

Next,  the  ratio  of  the  global  mole  numbers  of 
carbon  to  diatomic  hydrogen  is  calculated. 


f  .  Vf  / 

C  VC  '  VH, 


(3-2-2) 


This  ratio  is  also  constant  during  the  coabustion 


process 


An  initial  estimate  of  the  mole  fraction  of 


diatoaic  oxygen  is  required  and  is  generally 
initialized  as, 

EST(  o) 

°2 

At  this  point,  the  nole  fraction  of  atonic  oxygen  is 
calculated, 

%  -  K°  iRf  (3.2-3) 

followed  by, 

A  ■  (K°HVn0^T  ‘,1  KH)*(l  ♦  fc/2^  (3-2-4) 

»  -  .  >)•(*  *  ’c)  0-2*5) 

Ihe  sole  fraction  of  diatoaic  hydrogen  is  new 
determined. 


which  allows  the  aole  fraction  of  water  to  fellow, 

nH20  “  KH2°  n»2  iRf  (312-7) 

and  the  sole  fraction  of  OH, 


r* 

i 

<*“* 

32 

o 

as  well  as  ator.ic  hydrogen. 

*.•> 

n  -  kk  Jnj" 

H  *  H2 

(3-2-9) 

Mow  the  global  sole  fraction  of  diatomic  hydrogen 

is  calculated. 

— 

^H*  ■  nH2  4  nH20  *  (n0H  * 

(3-2-10) 

followed  by  the  global  sole  fraction  of 

carton. 

f 

„c  _  s 

nc  -  nK2  •  fC 

(3-2-11) 

*  •> 

Next,  the  sole  fractions  of  CO  and 

CO  are 

* 

■ 

obtained. 

1 1 

"CO  ■  "c  /  (>  ♦  *C°2  V%f  ) 

(3-2-12) 

*  \ 

s 

• 

nC02  ■  nc  -  nco 

(3-2-13) 

1  * 

and  finally  the  global  sole  fraction  of 

diatonic 

'•  y- 

oxygen , 

g  EST  /  \  , 

(3-2-lh) 

n°2  "  n°2  *  nC02  00  *  1,0  *  n°H  *  nH2°  >  2 

Using  the  above  sole  fractions  that 

are  based  on 

rS 

the  estimated  sole  fraction  of  diatonic 

oxygen,  a 

; 

calculated  value  of  *«i  is  deternined. 

calc  g  ,  g 

t 

f0  “  n0  /nH 

°2  °2  H2 

(3-2-15) 

J 

: 

i  n 

*  *-> 

• 

s 

R  C 

t  -• 

If  the  absolute  value  of  *o2*rox  is  large 
(greater  than  .0000001  for  exaaple)  a  new  estiaate  of 
the  aole  fraction  of  diatonic  oxygen  is  Bade  using  the 
following  empirical  relation. 


EST(n+l)  EST(n)  ,  calc  v* 

>o2  ■  «o2  •  (V  S  ) 


(3-2-16) 


and  the  iteration  is  started  again  by  returning  to 
equation  (3.2-3). 

For  gas  aixtures  involving  only  hydrogen  and 
oxygen,  the  above  iteration  can  be  used  by  specif  ying 
that  *c  is  equal  to  zero. 

For  gas  aixtures  containing  both  carbon  and 


and  is  inserted  after  equation  (3.2-21) .  A  calculated 


value  of 


V°a 


is  deteraiaed  as. 


calc 

rH2/02 


*  /  * 
**2  /  »»02 


(3-2-23) 


which  is  inserted  after  equation  (3.2-15). 

If  the  absolute  value  of  r„^-  *£0,  large 


(greater  thaa  .0000001  for  example)  then  the  estimate 
of  the  sole  fraction  of  diatoaic  nitrogen  is  improved 


by  the  following  espirical  relation. 


EST(n+l) 


EST(n) 


calc  •* 

/  V°2  ) 


(3-2-24) 


This  relation  is  placed  after  equation  (3.2-16)  which 
improves  the  previous  estimate  of  the  sole  fraction  of 
diatoaic  oxygen.  If  the  gas  aixture  contains  only 
hydrogen,  oxygen  and  nitregen,  the  value  of  fcis  set 
egual  to  zero. 

*htn  the  diluent  component  of  the  gas  aixture  is 
an  inert  specief^the  composition  iteration  is  started 
with  deteraination  of  the  ratio  of  the  global  mole 
number  of  the  inert  gas  to  that  of  the  diatomic 


hydrogen. 


*  *  /  8 

fx  "  v*  /  VH, 


(3-2-25) 


This  is  in  turn  followed  by  calculation  of  r0  from 

2 

eguation  (3.2-1)  and  an  initial  estimate  of  the  mole 
fraction  of  diatoaic  oxygen, 

EST(o) 

no2  ”  -1 

The  aole  fraction  of  atomic  oxygen  is  deterained 
froa  egaution  (3.2-3).  This  is  in  turn  followed  ty 


xv  7^  v*  m  "JV  :.%■  :v.  "vr  v  v-  * 
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calculation  of  the  following  two  quantities, 

A  «  (K°HV^  ♦  KH  )  •  (l  ♦  (3-2-26) 


B  -  (KH2° 


1  )  •  (l  ♦  fx  )  (3-2-27) 


The  sole  fractions  of  atoaic  hydrogen,  water  and  OH  are 
obtained  fcon  equations  (3.2-7)  to  (3.2-9)  .  The 
global  sole  fraction  of  diatonic  hydrogen  is  obtained 
fron  eg.  (3.2-10).  The  global  aole  fraction  of 
diatonic,  oxygen  is  obtained  fron. 


EST 


(""i 


2  u2  \  "2V 

•ade 


-  +  %)/2 


(3-2-28) 


,0  ^  ''OH  T  ''Qj 

The  value  of  is  determined  fron  equation 

(3.2-15).  A  new  estimate  of  the  nole  fraction  of 

diatonic  oxygen  is  aade  fron  eguation  (3.2-16).  The 

calc 

iteration  is  continued  until  fo  “  is  small. 

2  2 

The  procedure  described  in  section  3.1  for 
calculation  of  the  detonation  parameters  of  a  stable 
detonation  wave  was  used  to  determine  the  parameters  of 
the  following  gas  mixtures, 

*o2  +  h2  +*co2 

*02  +  H2  +  N2 


i02  +  H2  +  He 


i02  ♦  H2  ♦  Ar 


& 


*  _  •  m  •  •  • 

•  „  •  ,  »  „  • 
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*"/*r  ,*•*•*' 
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v-vy; 
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The  detonation  paraaetets  of  *  Wj  and  the 

gas  speed  were  calculated  for  different  initial 
teaperatures  of  100,  20C,  300,  400  and  500  degrees 
Kelvin.  Xhe  initial  pressures  considered  were  .1,  .5, 
1.0,  2.0,  and  5.0  ataospheies.  Ihe  results  of  these 
calculations  are  coapiled  in  Tables  (3-10)  to  (3-13). 

3.3  Initial  Pressure  7s.  Initial  Temperature. 

The  results  obtained  by  calculating  the  paraaeters 

of  a  stable  detonation  wave  describe  only  the  steady 

state  notion  of  the  wave.  It  is  therefore  difficult  to 

apply  then  in  explicitly  calculating  the  induction 

distance.  However,  it  is  interesting  to  notice  that 

P3 

the  pressure  ratio,  _  ,  nay  be  related  to  the 

pi 

magnitude  of  this  distance.  By  observing  plots 
(Figures  (3-1)  to  (3-4)  for  exanple)  of  this  ratio  as  a 
function  of  the  initial  gas  aixture  teaperature,  , 

we  see  a  relatively  sharp  increase  as  the  initial 
teaperature  decreases. 

Fcon  the  energy  equation  considered  before  and 
after  the  detonation  wave,  it  can  be  shown  that  the 
diaensionless  relative  energy  transfer  froa  the  burned 
gas  behind  the  detonation  wave  to  the  shocked,  tut 


unburnt  gas  ia  the  region  between  the  shock  wave  and 
the  coabustion  wave  is. 


Ah*  /p,  \ 

Vi  "  vpj' 


(3-3-D 


A  complete  derivation  of  equation  (3.3-1)  is  given  in 
the  appendix. 

Equation  (3.3-1)  shows  that  the  diaensionless 

relative  energy  transfer  is  directly  proportional  to 

p3 

the  pressure  ratio  j-  ,  and  inversely  proportional  to 

the  initial  temperature  of  the  gas  mixture.  This 

result  is  reasonable  since  a  combustible  gas  mixture 

will  generally  have  a  higher  enthalpy  at  a  higher 

temperature.  This  in  turn  lowers  the  amount  of 

relative  energy  transfered  prior  to  tne  formation  of 

the  stable  detonation  wave.  Furthermore,  the 

definition  of  the  detonation  induction  distance  is 

relative  to  whether  we  are  considering  the  point  where 

the  uustable  wave  forms,  cr  the  point  where  it  has 

become  stable.  Therefore,  we  can  justifiably  relate 

the  relative  energy  transfer  to  this  distance. 

It  is  particularly  interesting  to  note  that  the 
P3 

ratio  —  ,  does  s.  >t  change  substantially  as  the 

pi 

initial  pressure, ,  is  changed.  This  result  allows 
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us  to  conclude  that  initial  temperature,  rather  than 
initial  pressure  has  a  aore  important  effect  on  the 
length  of  the  induction  distance.  However,  the 
induction  distance  does  change  markedly  with  changes,  in 
initial  pressure. 

It  is  interesting  to  observe  that  by  aultiplying 
both  sides  of  eguation  (3.3-1)  by  R1t1  we  obtain. 


‘  **  (4  -  *) 


(3-3-2) 


Eguation  (3.3-2)  is  the  relative  energy  transfer 
per  unit  mass  of  gas  aixture.  Of  course,  a  plot  of 
this  parameter  (for  different  initial  pressures)  as  a 
function  of  the  initial  teaperature  would  show  tne  same 
trend  as  that  in  Figures  (3-1)  to  (3-4)  and  would  allow 
us  to  draw  the  sane  conclusions. 

However,  by  aultiplying  both  sides  of  eg.  (3.3-2) 
by  the  initial  density  the  relative  energy  transfer  per 


unit  voluae  is  obtained 


By  observing  values  of  this  parameter  (Table  (3-14)) 
for  the  hydrogen/oxygen/nitrogen  gas  aixture  for 
exaaple,  this  parameter  is  observed  to  increase 
dramatically  both  as  initial  pressure  is  increased  aai 
initial  tesperature  decreased.  Because  the  relative 
energy  transfer  per  unit  voluae  is  inversely 
proportional  to  initial  teaperature  and  proportional  to 
initial  pressure,  it  is  therefore  proportional  to  the 
initial  density  of  the  gas  aixture. 
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CHAP1EB  FCUfi 


EXPEBIWENIA1  ANALYSIS 


In  order  to  deternine  the  effect  of  the  speed  of 
sound  and  the  initial  density  on  the  induction 
distance,  several  hydrogen/oxygen/dilueut  gas  aixtures 
initially  at  tooa  temperature  (approxiaately  300 
degrees  K)  and  atmospheric  pressure  were  ignited  in  a 
o.4  aeter  long  combustion  tube. 

Ihe  proportion  and  types  of  diluent  species  used 
were  nitrogen  (  n2  ) ,  argon  (Ac),  heliua  (He)  and  carbon 
dioxide  (  co2  ).  The  aole  nuabers  cf  the  hydrogen  and 
oxygen  were  one  and  one  half  respectively.  Originally, 
the  aole  number  for  the  carbon  dioxide  was  one,  but  it 
was  lowered  to  one  half  after  it  was  discovered  that 
its  detonation  induction  distance  was  longer  than  the 
length  of  the  coabustion  tube  for  the  given  initial 
conditions. 

Also,  the  hydrogen/oxygen/heliun  and 
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hydrogen/oxygen/argon  gas  aixtures  were  burned  at 
atnospheric  pressure  and  an  initial  tenperature  of  140 
degrees  Kelvin. 


4.1  Apparatus 


The  cylindrical  combustion  tube  in  the  experiment 
had  au  inside  diaaeter  of  five  centimeters.  Along  it 
(see  figure  (4*1)}  was  an  array  of  thirteen  randomly 
spaced  holes  into  which  mounts  containing  piezoelectric 
quartz  pressure  transducers  could  be  placed.  The 
entire  tube  was  encased  in  a  steel  cooling  jacket 
through  which  liquid  nitrogen  could  be  passed  in  order 
to  cool  the  gas  aiztures  before  ignition  in  the  low 
temperature  part  of  this  experiment. 

The  gases  for  the  experiment  were  obtained  from 
commercially  available  cylinders  and  were  of  industrial 
purity.  The  individual  gases  were  premixed  in  a  remote 
mixing  chamber  before  iujection  into  the  combustion 
tube.  Injection  of  the  gases  was  accouplished  tnrough 
an  orifice  in  the  ignitor  assembly  (see  Figure  (4-2)). 
The  mixture  ratios  were  controlled  through  a  series  of 
regulators,  valves  and  gauges  contained  within  a 
control  panel  (see  figure  (4-3) } . 

After  injection  of  the  gas  mixture  into  the  tube, 
both  ends  were  closed  by  use  of  valves.  The  mixture 
was  then  ignited  at  one  end  by  passing  an  electrical 
current  through  a  thin  piece  of  nichrome  wire  wrapped 


across  two  electrodes  in  the  ignitor. 


4.2  Data  Seduction. 


Have  speed  was  aeasured  at  each  of  the  thirteen 
locations  along  the  tube  using  the  pressure  transducers 
in  conjunction  with  a  dual  bean  oscilloscope.  Two 
successive  transducer  locations  were  monitored  per  run. 
The  transducer  closest  to  the  ignitor  was  used  to 
trigger  the  pair  of  beaas  on  the  oscilloscope.  The 
wave  speed  was  deterained  by  first  oeasuring  the 
distance  that  the  second  (lower)  bean  travelled  before 
the  wave  passed  the  second  transducer.  Passage  of  the 
wave  was  accompanied  by  an  abrupt  junp  in  the  trace  due 
to  the  aoaentary  pressure  increase.  This  distance 
corresponded  to  the  wave  tine  between  traoducers.  inis 
tiae  was  taen  divided  into  the  distance  betweea  the 
transducers  which  yielded  the  average  wave  speed. 

Because  the  wave  passage  was  on  the  order  of 
microseconds,  a  photograph  of  the  oscilloscope  beaus 
(see  Figure  (4-4)  for  example)  traversing  the  screen 
was  required  for  each  run.  Ihe  tvo  transducers  were 
successively  repositioned  until  a  complete  profile  of 
the  average  wave  velocity  along  the  tube  had  been 


obtained. 

4.3  Experimental  Results. 

The  detonation  induction  distances  for  each  of  the 
gas  aixtures  were  determined  by  observing  the 
approximate  (plus  or  ainus  IS  cn.)  position  cn  the 
velocity  profiles  where  the  overshoot  described  in 
section  1.1  occurred.  Have  velocities  obtained  are 
coapiled  iu  Tables  (4-1)  to  (4-6)  and  these  values  are 
plotted  in  Figures  (4-6)  to  (4-11). 

It  is  interesting  to  note  that  the  induction 
distance  for  the  hydroges/oxygen/nitroyen  and  the 
hydrogen/oxygen/heliua  gas  aixtures  are  identical  for 
the  initial  conditions  used  in  tbrs  experiment.  This 
observation  is  consistent  with  the  data  of  Bellinger, 
Fong  and  Else  (ref.  28)  who  found  that  for  the  same 
initial  conditions  but  only  half  the  aaount  of  dilueut, 
the  induction  distances  were  the  same.  The  distances 
obtained  by  Bollinger,  Fong  and  Edse  are  less  than 
those  found  in  this  experiaent  which  is  to  be  expected 
since,  in  general,  the  greater  the  amount  of  diluent, 
the  longer  the  induction  distance. 

To  calculate  the  initial  Each  nuaber  of  the  flame 
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front,  flaae  speeds  obtained  by  Jones  (ref.  31)  for 
each  of  me  Mixtures  were  used.  The  Macb  nuatec 
relationship  is, 

h 

Mf  =  uf/(  y1  Rx  Tx)  (4. 3-1) 

The  induction  distances,  initial  flaae  front  Mach 
numbers  and  flame  speeds  are  compiled  in  Table  {4-7) 
where  the  data  of  Bollinger,  fong  and  £dse  have  also 
been  included.  Table  (4-8)  contains  the  data  from  the 
experiments  at  low  initial  temperature. 


CHAPTEH  FIVE 


CONCLUSIONS 


Sven  though  there  is  uuch  evidence  available  that 
confines  the  effect  of  both  initial  temperature  an-d 
initial  pressure  on  the  induction  distance  of  a 
combustible  gas  mixture,  it  is  inherently  difficult  to 
draw  conclusions  as  to  the  relative  importance  between 
these  two  parameters.  This  is  particularly  true  if 
they  are  to  be  based  solely  on  experimental  data. 
Therefore,  it  would  be  convenient  to  determine  a 
theoretical  relationship  that  would  clearly  show  which 
is  the  sore  important  factor. 

In  this  study,  the  dimensionless  relative  energy 
transfer  to  the  gases  behind  the  head  of  a  statla 
detonation  wave  has  been  analyzed  theoretically  in 
detail  for  several  hydrogen/oxygen/diluent  gas 
mixtures.  Plots  of  this  quantity  {Figures  (3-1)  to  (3 
4))  as  a  function  of  the  initial  temperature  for 
various  initial  pressures  show  that  the  initial 
temperature  appears  to  be  the  more  significant  factor 
in  changing  this  parameter. 
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In  fact,  the  aaxiaua  change  in  the  relative  energy 
transfer  is  only  on  the  order  of  ten  percent  when  the 
initial  pressure  is  increased  by  a  factor  of  fifty. 

The  aaxiaua  difference  occurs  when  the  initial 
teaperature  is  very  low  <100  degrees  Kelvin  in  this 
study)  and  decreases  rapidly  as  the  initial  teaperature 
is  increased.  These  results  are  based  on  the  gas 
aixtures  considered  in  this  study  but  there  is  no 
reason  to  believe  that  other  coabustible  gas  aixtures 
would  not  behave  in  this  Banner.  Considering  the 
relative  energy  transfer  per  unit  aass  (equation 
(3.3-2))  results  in  the  saae  conclusion. 

Shea  the  relative  energy  transfer  per  unit  voluae 
is  plotted  as  a  function  of  the  initial  teaperature  for 
various  initial  pressures  (figure  (3-5)  for  exaaple) , 
this  paraaeter  is  observed  to  change  markedly  with 
initial  pressure  as  well  as  with  initial  teaperature. 
Therefore,  no  insight  into  the  relative  importance 
between  teaperature  and  pressure  is  obtained, 
particularly  in  describing  detonations  in  fixed 
voluaes.  This  is  because  in  a  fixed  voluae,  changes  in 
density  neccessitate  changes  in  aass.  lhis  result 
leads  back  to  consideration  of  the  diaensionless 
relative  energy  transfer  or  the  relative  energy 


-J 
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transfer  per  unit  aass. 

As  far  as  applying  the  relative  energy  transfer  to 
the  detonation  induction  distance,  because  it  is 
proportional  to  initial  temperature,  it  is  easy  to  say 
that  it  is  therefore  proportional  to  the  .induction 
distance.  However,  there  are  two  major  disagreeaents 
with  this  conclusion.  The  first  arises  because  the 
induction  distance  for  two  different  gas  aixtures  can 
be  the  sane  but  the  relative  energy  transfer  nay  not. 
The  hydrogen/oxygen/helium  and  hydrogen/oxygen/nitrogen 
systems  are  an  example.  The  second  conflict  arises 
when  two  different  gas  aixtures  with  the  same  relative 
energy  transfer  do  not  have  the  saae  induction 
distance.  The  hydrogen/oxygen/helium  and 
hydrogen/oxygen/argon  systems  are  examples  of  this 
case. 

It  is  therefore  difficult  to  say  whether  a 
relationship  for  the  length  of  the  induction  distance 
in  a  combustible  gas  aixture  can  be  obtained  based  on 
the  concept  of  relative  energy  trausfer  to  the  staLle 
detonation  wave. 

Hhen  the  initial  flaae  front  Sach  nuaber  and  its 
relationship  with  the  induction  distance  is  studied, 
one  major  problem  arises.  The  relationship  is 


complicated  by  the  fact  that  different  gas  mixtures  can 
have  the  same  induction  distance  but  different  initial 
flase  front  Bach  nuabers.  From  the  data  obtained  in 
this  study  and  the  data  compiled  by  other 
investigators,  it  is  seen  that  gas  sutures  with 
relatively  high  flase  front  Bach  nuabers  can  have 
either  long  or  short  induction  distances.  The 
hydrogen/oxygen/argon  and  hydrogen/oxygen/carbou 
dioxide  aixtures  display  this  effect. 

Therefore,  consideration  of  the  relative  energy 
transfer  and  the  speed  of  sound  in  the  unburnt  gas 
mixture  show  these  parameters  to  be  inadequate  in 
themselves  to  describe  the  behaviour  of  the  detonation 
induction  distance  in  cosbustible  gas  aixtures. 

The  combustion  tube  used  in  this  study  did  not 
have  transducer  mounts  in  the  region  very  close  to  the 
ignitor.  This  presents  a  problem  in  judging  the 
accuracy  of  the  value  obtained  for  the  induction 
distance  in  each  of  the  low  temperature  hydrogen/ 
oxygen/inert  gas  aixtures.  This  is  particularly  true 
in  the  case  vhere  the  inert  species  vas  argon. 

it  must  be  mentioned  that  the  wave  speeds  after 
detonation  for  the  low  temperature  mixtures  do  not 
agree  veil  vith  the  calculated  values  for  the  stable 


detonation  waves.  This  is  observed  by  looting  at  the 
profiles  of  Figures  (4 -1U)  and  (4-11)  and  cospacing  the 
wave  speeds  after  detonation  with  these  coapiled  is 
Tables  (1-12)  and  (3-13).  The  calculated  values  are 
consistently  higher  than  what  was  actually  seasured. 

To  iaprove  the  experiaental  results  of  this  study, 
it  is  suggested  that  a  coabustion  tune  of  appcoxiaately 
ten  feet  in  length  he  built  containing  tcansducer 
locations  at  three  inch  intervals  froa  the  ignitor.  In 
additica  to  construction  of  a  new  tube,  the  use  of  at 
least  four  transducers  pec  run  is  advised  such  that 
aore  data  can  be  obtained  in  a  shorter  aaount  cf  tine. 
In  general,  the  use  of  aany  transducer  acunts  along  the 
length  of  the  coabustion  tube  will  produce  aucii  aore 
coaplete  profiles  and  allow  aore  accurate  estiaation  of 
the  detonation  induction  distance. 


TABLE  3-1 


Equations  Foe  Eguilitriua  Constants 

K1 

=2 . 42- T'125.  «p  {(^0  -  A  } 

8°-%o,T'*  *  “■“{'&-)  A|  ’p* 


“  *0H  *  exp 


K^")oh  Al 


0  '  So  •  •*»  H¥)  a} 


TABLE  3-2 


Coefficients  Of  Equilibrium  Constants 


4*2**  X 

0 

*°2*H2'-H2° 

16.351 

41.289 

5.8457 

.08531 

16.405 

41.20? 

5.8133 

.08483 

16.451 

41.  112 

5.7810 

.08440 

16.495 

40.993 

5.7490 

.08403 

16.543 

40.862 

5.7171 

.C8370 

16.585 

40.719 

5.6854 

.08352 

16.602 

40.553 

5.6594 

.08328 

16.620 

40.391 

5.6334 

.08307 

16.635 

40.225 

5.6076 

.08295 

16.624 

40.061 

5.5820 

.08283 

16.645 

39.901 

5.5564 

.03284 

16.645 

39.738 

5.5310 

.03232 

16.638 

39.575 

5.5056 

.08283 

16.625 

39.413 

5.4804 

.06286 

16.612 

39.256 

5.4553 

.03929 

16.602 

39.087 

5.4303 

.08300 

table  3-3 


4000 


Coefficients  Of 

Equilibrium 

A 

Constants 

**2—  N 

tOj*CO“C02 

39.865 

4.4620 

1.5090  -5 

39.865 

4.4651 

1.5137  -5 

39.861 

4.4669 

1.5199  -5 

39.845 

4.4687 

1.5236  -5 

39.809 

4.4696 

1.5300  -5 

39.734 

4.4/05 

1.5366  -5 

39.701 

4.4709 

1.5441  -5 

39.666 

4.4705 

1.5526  -5 

39.630 

4.4696 

1.5615  -5 

39.586 

4.4687 

1.5703  -5 

39.539 

4.4669 

1.5786  -5 

39.488 

4.4647 

1.5870  -5 

39.428 

4.4620 

1.5956  -5 

39.365 

4.4564 

1.6040  -5 

39.302 

4.4549 

1.6128  -5 

39.220 

4.4513 

1.6214  -5 

TABLE  3-4 


°K) 

00 

00 

00 

00 

00 

00 

00 


30 

00 

00 

00 

00 

00 

00 

00 


Seduced  Sensible  Enthalpies 


H2 

H 

°2 

0 

3.799 

2.5 

4.  188 

2.543 

3.819 

2.5 

4.208 

2.542 

3.839 

2.5 

4.227 

2.541 

3.  859 

2.5 

4.245 

2.540 

3.878 

2.5 

4.263 

2.539 

3.897 

2.5 

4.281 

2.538 

3.916 

2.5 

4.298 

2.537 

3.934 

2.5 

4.315 

2.537 

3.952 

2.5 

4.332 

2.537 

3.969 

2.5 

4.348 

2.53b 

3.  987 

2.5 

4.364 

2.536 

4.003 

2.5 

4.379 

2.536 

4.020 

2.5 

4.395 

2.537 

4.037 

2.5 

4.409 

2.537 

4.053 

2.5 

4.424 

2-537 

4.069 

2.5 

4.438 

2.538 

OH 

3.877 

3.894 

3.911 

3. 928 

3.944 

3.959 

3.974 

3.989 

4.003 

4.017 

4.031 

4.044 

4.056 

4.069 

4.081 

4.093 


TABLE  3-5 

Seduced  Sensible  Enthalpies 
T 

s  k  * 


(^T^i 


H20 

N2 

NO 

N 

CO 

C°2 

5.238 

3.992 

4.114 

2.500 

4.026 

6.316 

5.286 

4.008 

4.  128 

2.501 

.  4.042 

6.358 

5.333 

4.024 

4.141 

2.501 

4.057 

6.398 

5.378 

4.038 

4.153 

2.501 

4.071 

6.435 

5.421 

4.052 

4.165 

2.502 

4.0  85 

6.470 

5.463 

4.066 

4.176 

2.502 

4.098 

6.504 

5.503 

4.078 

4.187 

2.503 

4.110 

6.536 

5.542 

4.090 

4.  197 

2.504 

4.122 

6.566 

5.580 

4.102 

4.207 

2.505 

4.133 

6.595 

5.616 

4.113 

4.  2 16 

2.506 

4.144 

6.62J 

5.651 

4.  124 

4.226 

2.500 

4.  154 

6.649 

5.665 

4.134 

4.234 

2.509 

4.164 

6.674 

5.717 

4.  144 

4.24  3 

2.511 

4.  174 

6.690 

5.749 

4.154  . 

4.251 

2.513 

4.183 

t.  722 

5.779 

4.163 

4.258 

2.515 

4.192 

6.744 

5.809 

4.  172 

4.266 

2.518 

4.200 

6.7b5 

TABLE  3-7 


Diaensioaless  Specific  Heats 


T(°K) 

H2 

H 

2500 

4.315 

2.5 

2600 

4.347 

2.5 

2700 

4.378 

2.5 

2800 

4.407 

2.5 

2900 

4.433 

2.5 

3000 

4.458 

2.5 

3100 

4.484 

2-5 

3200 

4.510 

2.5 

3300 

4.535 

2.5 

3400 

4.560 

2.5 

3500 

4.564 

2.5 

3600 

4.609 

2.5 

3700 

4.632 

2.5 

3800 

4.656 

2.5 

3900 

4.679 

2.5 

4000 

4.701 

2.5 

4.681 

2.508 

4.707 

2.509 

4.733 

2.511 

4.758 

2.513 

4.782 

2.516 

4.806 

2.518 

4.829 

2.521 

4.851 

2.525 

4.872 

2.529 

4.893 

2.533 

4.913 

2.537 

4.931 

2.541 

4.949 

2.546 

4.966 

2.551 

4.982 

2.557 

4.998 

2.562 

OH 

4.339 
4.333 
4.  360 
4.381 
4..  400 
4.418 
4.435 
4.452 
4.467 
4.481 
4.495 
4.508 
4.521 
4.533 
4.545 
4.556 


TABLE  3-10 


Detonation  Parameters  For 


io2  *  h2  ♦  ico2 

At  . 1  And  5  Atmospheres  Initial  Pressure 
—  *1  atm • 


r1(°K) 

t3(°k) 

P3(atm) 

W^m/s) 

CO 

> 

<*> 

( m/s ] 

100 

3250.52 

4.  28 

1926.87 

1026.53 

900.34 

200 

3168.72 

2.  10 

1897.46 

1021.87 

875.59 

300 

3127.24 

1.39 

1875.65 

1021.25 

854.40 

400 

3101.54 

1.03 

1856.93 

1023.44 

833.49 

3082-58 

-82 

1838.09 

1025.80 

812.29 

p i  9  5  atm. 

100 

3813.00 

231.8 

2014.50 

1079.07 

935.43 

200 

3733.38 

114.6 

1992.25 

1077.69 

914. 5o 

300 

3695-45 

76.0 

1975.94 

1080.03 

895-91 

3674.75 

56.9 

1962.31 

1083.60 

878.71 

500 

3659.86 

45.3  • 

1948.11 

1087.67 

660.44 

1A8LE  3-11 


Detonation  Parameters  For 


i02  ♦  H2 

+  N2 

1  Hud  5 

Atnospheres  Initial 

Pressure 

P1  «  .1  atm. 

t3(°k) 

P3(atm) 

w1(m/s) 

( m/s ) 

U3(m/s) 

3076.23 

4.92 

2167.68 

1162.83 

100U.9 

3015.84 

2-41 

2134.08 

1 15a. 05 

980.0 

2986.41 

1.59 

2109.01 

1150.98 

958.0 

2968.67 

1.  18 

2085.90 

1150.08 

935.8 

2957.45 

.94 

2066.37 

1150.41 

915.9 

Px  *  5  atm. 

3499.02 

265.0 

2281.36 

1236.06 

1045.3 

3457.83 

133.0 

2255.86 

1230.62 

1025.2 

3442.43 

88.  1 

2236.92 

1229.8o 

100 7.  1 

3436.98 

65-8 

2220.38 

1230.96 

989.4 

3435.83 


52.4  2204.57  1232.85 


971.7 


TABLE  3-13 


Detonation  Parameters  for 
i02  +  Hg  ♦  Ar 

1  And  5  Atmospheres  Initial  Pressure 


r3(°K) 

Pi  *  .1 

P3(atm) 

atm. 

w^m/s) 

W3(m/s) 

U^m/s) 

3253.20 

5.35 

2019.8  0 

1087.54 

932.26 

3165.59 

2.61 

1981.90 

1074.61 

907.30 

3120.09 

1.71 

1954.70 

1068.76 

885-95 

3090.99 

1.27 

1931.61 

1065.70 

865.91 

3070.  14 

1.00 

1910-19 

1063.99 

84G.20 

3866.65 

p i  9  5  atm. 

303.6  2171.17 

1 184.70 

98b. 47 

3776. 10 

148.8 

2137.38 

1172.37 

965.01 

3730.89 

98.0 

2113.89 

1 1o  7.  14 

946.74 

3703.49 

72.9 

2094.30 

1164.73 

929.57 

3684.85 

57.8 

2076.39 

1 163.68 

912.71 

TABLE  4-1 


Have  Velocities,  ,  Vs.  Distance  Fcoe  Ignitor 

io2  +  h2  ♦  ico2 

(1  Ate.  Initial  Pressure  and  °300  K  Initial  Temp.) 


Distance  (a) 

Han  *1  (e/s) 

Bun  1 2  (a/s) 

Bun  *3  (e /£) 

1.49 

162 

162 

162 

2.41 

941 

941 

941 

3.01 

1360 

1360 

1580 

3-31  • 

1613 

1741 

16  13 

3.94 

2511 

2651 

251 1 

4.53 

2162 

2179 

2214 

•  4.84 

2191 

2125 

2093 

5-77 

2169 

2126 

2105 

6.23 


2044 


2007 


2032 


TABLE  4-2 

Have  Velocities,  wi  ,  Vs.  Distance  Froa  Ignitor 
*o2  +  H2  +  N2 

(1  Atn.  Initial  Pressure  and  °300  K  Initial  leap.) 
Distance  (a)  Bun  *1  (a/s)  Bun  #2  (o/s)  Bun  *3  (a/s) 


TABLE  4-3 


Save  Velocities, Wj  ,  vs.  Distance  Pros  Ignitor 
*02  ♦  H2  ♦  He 

(1  Ate.  Initial  Pressure  and  °300  K  Initial  leap.) 
Distance  (a)  Run  *1  (a/s)  Bun  12  (a/s)  Bun  *3  in/s) 


TABi.£  4-4 


Nave  Velocities,  ,  Vs.  Distance  From  Ignitor 

i02  +  H2  +  Ar 


(1  At  B. 

Initial 

Pressure 

and  °300  K  Initial 

leap.) 

Distance 

(a)  Ban 

*1  (a/s) 

Bun  *2  (b/s)  Bun 

#3  (B/S) 

1.49 

312 

298 

319 

1.80 

960 

1131 

1002 

2.41 

2674 

2508 

2362 

3.01 

2347 

240b 

23t>7 

3.31 

2336 

2336 

2336 

3.94 

2272 

2272 

2236 

4.53 

2179 

2197 

2179 

4.84 

2262 

2262 

2262 

5.77 

2236 

2147 

2147 

TABLE  4-5 


Have  Velocities,  wj_.  Vs.  Distance  Froa  ignitor 
i02  +  H2  ♦  He 


(1  At  a. 

Initial 

Pressure 

and  0 140  K  Initial 

leap.) 

Distance 

(a)  Bun 

#1  (a/s) 

Bun  *2  i a/s)  Bun 

#3  ia/s) 

.43 

280 

286 

320 

.90 

570 

545 

600 

1.49 

868 

995 

1194 

1.80 

76  2 

469 

554 

2.41 

1051 

— 

— 

3.01 

2172 

2389 

2275 

3.31 

1353 

1449 

966 

3.94 

1296 

— 

— 

4.53 

1003 

841 

1152 

* 


Data  not  usable 


TABLE  4-6 


Have  Velocities,  ,  vs.  Distance  Pros  Ignitor 

J02  +  H2  ♦  Ar 

(1  At«.  initial  Pressure  and  °140  K  Initial  leap.) 


Distance  (a) 
.43 
.90 
1-49 
1.80 
2.41 
3.01 
3.31 


Sun  #1  (a/s) 
609 
2213 
1843 
442 
2121 
2120 
2018 


£un  #2  (a/s) 
563 
2326 
1393 
1637 
1840 
2222 
1945 


fiun  #3  (a/s) 

1528 

1707 

2015 


* 


Data  not  usatle 


TABLE  4-7 


Flaae  Speed,  Speed  of  Sound,  Initial  Flaae  Front 

Hach  Nuaber,  Initial  Density  And  Induction  Distance 
approx.  300  °K  and  p1  *  1  atm. 


Diluent 

uf  (»/s) 

Wa^«/S) 

M- 

f 

P1  (Kg/a**  3) 

(CB. 

N2 

4.60 

43  5.  62 

.0106 

.7479 

325 

Ar 

4.43 

398.68 

.0111 

.9418 

225 

He 

6.44 

646.87 

.0099 

.3578 

325 

*C°2 

6.07 

367.94 

.0147 

.8129 

375 

iN2 

7.89 

467.18 

.0168 

-6504 

228* 

^Ar 

8.00 

435.73 

.0184 

.7716 

179* 

Ar 

4.43 

398.68 

.0111 

.9418 

240* 

iHe 

9.34 

600.28 

.0156 

.4065 

228* 

*co2 

6.07 

412.54 

.0147 

.8129 

— 

*  Data  from  reference  28  (15  as  diameter  tube) 


TA8LE  4-8 


Flame  Speed,  Speed  o f  Sound,  Initial  Flaae  Front 

flach  Number,  initial  Density  And  Induction  Distance 

approx.  140  °K  and  pj  8  1  atm. 

Diluent  u^(h>/s)  wa^(a/s)  Pj(Kg/B**3)  d^(cm 

He  3.076  454.07  .0068  .7666  275 

Ar  2.117  275.85  .0076  2.018  75 
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I.  INTRODUCTION 

This  analysis  investigated  the  effect  of  heat  addition 
(due  to  fuel  mass  addition)  on  the  thermodynamic  and  gas- 
dynamic  properties  of  an  established  one-dimensional, 
chemically  reacting,  subsonic  flow.  In  particular  the  effect 
of  the  fuel  mass  addition  on  the  upstream  conditions  of  the 
primary  flow  (air  in  this  study)  was  examined. 

It  was  assumed  that  the  flow  was  adiabatic  relative  to 
the  environment  and  friction  was  negligible  between  the  flow 
and  the  walls  of  the  duct.  The  length  of  the  duct  from  the 
measuring  point  of  the  upstream  primary  flow  conditions  to 
the  point  where  complete  combustion  has  occured  (chemical 
equilibrium  prevails)  is  assumed  to  be  long  enough  to  allow 
for' secondary  flow  injections  (fuels),  and  complete  homo¬ 
genous  mixing  of  all  gases.  Most  importantly,  it  was 
assumed  that  at  some  point  past  the  combustion  zone,  the 
flow  becomes  choked.  This  last  assumption  allows  the  problem 
to  be  solved  and  in  practice  is  found  in  turbojet  operations. 

Only  the  theor  tical  analysis  is  included  in  this  report. 
Gaseous  hydrogen  has  been  used  as  the  fuel  because  of  its 
high  heat  release  and  availability.  Due  to  the  complications 
involved  in  building  an  experimental  device  which  closely 
approximates  a  basic  gas  generator,  the  experimental  results 
have  yet  to  be  completed. 


II.  THEORETICAL  ANALYSIS 
A.  General 

The  objective  of  this  report  was  to  theoretically 
simulate  the  effect  of  fuel  flow  changes  on  the  temper¬ 
ature,  pressure,  and  speed  of  the  flow  of  a  basic  gas  gener¬ 
ator  which  is  the  heart  of  the  turbojet,  turbofan  and 
turbo-prop  engine.  This  problem  has  been  approached  in 
past  literature  under  two  areas,  compressor-turbine  matching 
and  acceleration,  deceleration  problems  of  turbojets.  The 
latter  being  more  in  tune  to  this  report.  However  both  areas 
have  only  been  lightly  discussed  in  compressor  stall  mater¬ 
ial,  and  quantitative  information  on  them  is  scarce. 

Industry  has  concentrated  little  effort  towards  studying 
acceleration,  deceleration  stalls  from  fuel  flow  changes  for 
they  feel  that  the  acceleration  and  deceleration  schedules 
built  into  the  fuel  controls  of  their  production  engines 
sufficiently  handle  the  problem.  However  it  is  felt  by 
these  authors  that  pressure  disturbances  (from  an  increased 
fuel  flow)  are  generated  much  faster  than  these  schedules 
can  react.  Therefore  engine  failure  results  despite  having 
been  thought  to  be  free  of  such  a  failure.  As  will  be  shown, 
this  relatively  simple  assumption  gives  a  feasable,  realistic 
explanation  to  unexplained  engine  failures  experienced  by 
some  of  the  high  performance  turbojets  on  fighter  aircraft. 
The  problem  is  more  prevelent  today  since  we  require  engines 
to  operate  near  their  peak  performances. 
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B.  Compressor  Maps 


Examining  compressor  maps  of  an  engine,  one  can  see 
how  fuel  flow  controls  or  schedules  do  decrease  the  chance  of 
stalls.  Figure  1  shows  a  typical  compressor  map  which  lends 
to  understanding  the  operating  limits  of  compressors.  The 
pressure  rise  through  the  compressor  (P3/P2)  is  plotted 
against  an  altitude  corrected  air  mass  flow  rate  (m  02/62). 

The  notation  for  figures  1  and  2  is  given  below. 

Table  1.  Compressor  map  notation 

6  -  P/Ps  P  -  altitude  pressure  Ps  -  standard  press. 

0  -  T/Ts  T  -  altitude  temp.  T  -  standard  temp, 

m  -  mass  flow  rate  s 

(2) -  subscript  for  compressor  inlet 

(3) -  subscript  for  compressor  exit 

(4) -  subscript  for  turbine  inlet 
N  -  compressor  RPM's 

Under  an  ideal  acceleration,  an  increase  in  fuel  would 

increase  the  combustion  temperature  which  in  turn  would  in- 

Vi 

crease  the  work  done  the  turbine,  increasing  the  comp¬ 
ressor  RPM.  The  schedules  now  built  into  production  engines 
maintain  a  path  such  as  path  Y.  Along  Y  the  temperature  ratio 
does  not  Increase  so  drastically  that  stall  conditions  are 
reached. 

However  this  is  not  a  complete  picture  of  the  gas  gener- 

khe. 

ator.  Nothing  in  the  compressor  maps  accounts  for .pressure 
being  sent  upstream  to  Ww. 

combustion  chamber  from  the  compressor.  These  pressure  dis¬ 
turbances,  as  shown  in  this  report,  disrupt  the  flow  from  the 
compressor.  This  can  lead  to  stall  in  the  compressor, 

3 


Thus  the  entire  com¬ 


pressor  may  be  disturbed. 

The  arrow  in  figure  1  shows  how  the  stall  or  surge 
line  decreases  as  altitude  is  increased.  Most  likely  the 
pressure  disturbance  from  a  fuel  flow  change  also  decreases 
the  pressure  ratio.  Therefore  if  path  Y  is  again  followed, 
it  is  possible  that  it  now  passes  through  the  stall  zone. 

As  stated  earlier  it  is  believed  by  these  authors  that 
the  actual  path  the  engine  takes  is  not  like  Y  but  more  like 
path  Z  in  figure  2,  Here  the  temperature  ratio  increases 
much  faster  than  the  turbine  and  compressor  speed* 

When  the  engine  is  operating  at  a  higher  CPR  and  higher  RPM, 
as  are  todays  engines  for  maximum  performance,  this  Increases 
the  chances  of  suchan  occurance  for  one  would  be  operating 
in  the  upper  right  hand  vicinity  of  figures  1  and  2  where 
there  is  less  room  for  nU/uvcutlov,  , 


Surge  Line 


Decel.  Line 


Previous  surge  line 


Compressor  map  variation 


Heat  Addition  Cnly 


r 

w 

Before  studying  flow  changes  due  to  fuel  mass  addi¬ 
tion  much  insight  into  the  problem  is  obtained  by  first 
looking  at  how  heat  addition  alone  causes  changes  to  an 
unheated,  subsonic,  inviscid  flow.  A  one-dimensional, 
constant  area  duct  is  used  and  the  ^,0^  is  assumed  to  be 
calorically  perfect.  Figures  3  and  4  illustrate  the  sub¬ 
script  notation  used  in  this  particular  analysis. 

At  some  point  along  the  duct,  heat  is  added  to  the 
flow  producing  a  local  change  in  pressure.  Since  the  flow 
is  subsonic,  the  pressure  disturbance  is  transmitted  both 
downstream  and  upstream  through  the  duct,  thus  altering 
the  initial  conditions  of  the  entire  flow.  In  particular 
the  pressure  waves  may  cause  the  static  pressure  of  the 
unheated  gas  to  increase  and  bring  about  a  reduction  ir.  the 
mass  flow  rate.  This  will  now  be  shown  to  be  true. 

Two  are  required.  The  first  is  that  the 

static  pressure  downstream,  that  of  the  heated  gas,  remains 

a/H/cC }  pressure. 

constant^therefore^  is  that  of  the  initial^  ?2  =  Pi).  The 
second  c<r> ictvtwm  is  that  the  compression  of  the  unheated 
gas  as  a  result  of  heat  addition,  is  isentropic  and  no 
energy  is  transferred  to  the  initial  gas  because  the  rate 
of  compression  is  moderate.  This  implies  that  =  p?  and 

0  o 

-l  3  Ti.  Thus  the  problem  has  only  3  unknowns,  M^,  Mg,  and 
Tg.  None  of  the  three  can  be  solved  for  directly.  Cne 


subscript  (1) 


► 


Pig.  3  -  Unheated  flow 


Pig.  4  -  Plow  after  heat  addition 


solves  the  problem  by  applying  the  continuity,  momentum, 
and  energy  equations. 


The  continuity  equation  written  for  figure  4  is 

PlUi  *  Pi Uj  - (1) 

Letting  pu  ■  pu  ■  ypM  ,  and  P2  =  Fj_ ,  equation  (1)  becomes 
**  ^RT 


FiKi  =  FjM2.  —(2) 

vr i  \fn 


With  no  work  being  performed  on  the  gas,  the  momentum 
equation  can  be  written 

Pi(1+YM?)  -  Pj/l+yMg).  - (3) 

The  energy  equation  for  this  flow  is 

O  O 

T2  =  Ti  +  a_ 

Cp 

or  Tj-T.d  ‘  — (*) 

where  T^  =  and  q  =  Joules  per  Kg  of  air,  Cp  *  j/Kg. K. 
So  q/CpTi  is  a  dimensionless  heat  release  factor.  The 
isentropic  relationships  give 


Px  »  Pi/l+lzi  M*Vl  - (5) 

T1  ■  Mi>>  T2  "  T2/(1+I=i  Ml>  -<6> 


Substition  of  equations  (4), (5),  and  (6)  into  (2)  gives 
a  relationship  between  two  unknowns  ,  y,2,  and  three 
k.nowns ,  Mi#  Tj,  and  q.  It  is  listed  below. 


A  second  relationship  between  and  Mg  i-s  obtained 
by  substition  of  equation  (5)  into  the  momentum  expression 
equation  (3).  This  gives  for  M£ 

M2  -  mN?-1  J  __(8) 


Equations  (7)  and  (8)  now  contain  only  two  unknowns 
and  Mg.  Neither  can  be  solved  for  explicity  so  some 
type  of  trial  and  error  iteration  must  be  used.  Cnee  Mi 
and  Mg  are  found  for  a  given  heat  input  (q/CpT°i),  the 
static  temperature  of  the  heated  gas  can  be  found  from 
equations  (4)  and  (6)  which  can  be  written  as 

T2  ’  Tl/*2\‘A+I^ 

w  v 


Also 


Y 

p2  *  Pjd+Y-1  M|)y-1 


Completing  the  analysis,  the  final  mass  flow  rate  is 
found  from  the  mass  flow  rate  general  expression. 


1 iZL  *  rr  m  y+1 
VYRT  (1+Y-l  MdilY-D^YRT0 


With  Tj  =  Tj_  and  F^  =  F]_  this  gives 


v*.  A  A.  A/.  A  Ay.  A  /.  AAA  AA^»  AA  A.  a  a 


These  expressions  show  that  heat  addition  to  a  flow  cus 
described  earlier,  increases  the  upstream  static  pressure^ 
slows  down  the  flow,  and  decreases  the  mass  flow  rate.  It 
is  easy  to  see  that  if  some  high  precision  device  such  as 
a  compressor  were  located  upstream  of  the  flow,  it  would 
react  to  the  heat  addition  downstream.  Figures  5-13  show 
the  effects  of  heat  addition  (q/CpT^)  to  an  inviscid  sub¬ 
sonic  flow  passing  through  a  constant  area  duct.  Equations 
(1)  -  (9)  were  used  to  obtain  these  figures  with  ?2  =  Pi  = 


8 


9 


10 


Mass  flow  rate  vs.  heat  addition 


Speed  of  heated  gas  vs.  heat  addition 


1 .  General 

The  effects  of  mass  addition  due  to  a  fuel  flow 
change  will  .now  be  discussed  in  this  section.  The  addi¬ 
tional  mass  in  the  flow  releases  heat  upon  combustion. 

Both  the  increased  heat  and  the  presence  of  the  additional 
mass  alter  the  flow  conditions  upstream  and  downstream. 

The  flow  -is  considered  to  be  steady,  one-dimensional, 
inviscid,  adiabatic  with  the  environment,  and  it  is  assumed 
that  complete  homogeneous  mixing  has  occured  before  combus¬ 
tion. 

Th*  jollo»i*3  anal  y/S  is  ba^eol  cr\  the  spec*'*/  condi'ti  oy\. 
that  at  some  point  past  the  combustion  zone  the  flow  Is 
choked.  This  is  a  realistic  assumption  for  turbine  nozzles 
are  often  choked  over  a  good  part  of  their  operating  range. 
Extending  this  analysis  one  can  see  that  with  a  choked 
afterburner,  the  turbine-afterburner  operation  of  the  turbo¬ 
jet  may  be  disrupted  by  way  of  fuel  changes  to  the  after¬ 
burner.  Therefore  had  choked  flow  not  been  assumed,  the 
problem  would  have  to  be  solved  with  incorrect  assumptions 
placed  on  the  combustion  gas.  {This  was  tried  and  produced 
inaccurate  results) 

Many  attempts  were  made  to  formulate  solutions  which 
accurately  simulated  fuel  flow  changes  and  their  effect  on 
the  flow.  Primarily,  modifications  to  the  previously  com¬ 
pleted  heat  addition  of  section  C  were  attempted.  But  these 


attempts  either  required  .non-reaiistic  assumptions  as 
mentioned  before,  cr  produced  equations  with  too  many  un¬ 
knowns.  Always  the  equations  of  motion  had  to  be  satisfied. 
The  difficulty  of  the  entire  problem  arises  because  the 
upstream  conditions  change.  One  can  not  simply  set  up  the 
equations  of  motion  by  adding  the  primary  and  secondary 
flows  then  solving  for  the  combustion  flow.  Because  the 
primary  flow  changes,  an  iteration  process  must  be  used. 

This  technique  is  outlined  in  the  next  subsection. 

2.  Solution  of  Problem 

The  addition  of  gases  to  a  primary  flow  through  a 
duct  may  be  accomplished  in  either  of  the  three  ways  shown 
in  figures  14  -  15.  The  secondary  flows  can  be  added  from 
the  side  and  may  or  may  not  increase  the  duct  area  (fig.  14) 
or  the  fuel  can  be  added  tangentially  (fig.  15).  the  setup 
in  fig.  15  was  used  for  it  simplified  the  momentum  equation 
and  in  the  future  would  make  the  physical  apparatus  easier 
to  construct.  It  is  realized  that  in  combustion  chambers 
the  mixing  of  the  fuel  and  combustion  occurs  almost  simul¬ 
taneously,  but  in  this  analysis  the  process  is  separated 
for  simplicity. 

Two  ways  can  be  used  to  get  the  values  of  the  combus- 

n'gorous  o»e 

tion  chamber  variables.  The  involves  using  equili¬ 

brium  constants  to  c«|eulc*fco  the  mole  fractions  of  each 
species  in  the  combustion  gas.  However  since  the  tempera¬ 
ture  in  turbojet  combustion  chambers  is  always  lower  than 
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2000°?',  dissociation  is  minimal  and  the  method  is  tedious. 

(For  a  complete  description  of  this  method  see  reference  3)« 
Instead,  dissociati'c  n  is  disregarded  and  the  co  w  has  £io  n 

is  assi/ned  fo  consist  oj  r\i  troge  n ,  n t  ajat^r 

t>£»  f>o  f  f  *a  »  d  cqrhon  dioxide  only. 

Both  methods  do  use  the  same  general  approach  in  cal¬ 
culating  the  speed  and  temperature  of  the  combustion  gas. 

This  involves  calculating  the  speed  from  two  different 
equations.  Cne  is  derived  from  the  momentum  and  continuity 
equations  and  the  other  derived  from  the  energy  equation. 

Both  require  an  estimated  combustion  temperature.  The  two 
speeds  are  compared  for  a  certain  temperature  and  a  new 
estimate  of  the  temperature  is  made  based  on  their  differ¬ 
ences.  This  iteration  is  continued  until  the  value  of  the 
speeds  from  each  equation  are  pruc  ti  c«  ''  yr  the  same.  This 
also  gives  the  correct  temperature.  Cnee  the  speed  and 
temperature  are  known  the  pressure  is  found  from  the  conti¬ 
nuity  equation. 

It  should  be  noted  at  this  point  that  problems  en¬ 
countered  most  often  when  programming  this  problem,  involved 
the  calculation  of  the  two  speeds  mentioned  above.  Figures 
16  and  17  give  qualitative  pictures  of  how  the  speeds  change 
with  temperature.  Referring  to  fig.  16,  for  an  initial 
temperature  estimate  points  A  are  found  from  each  equation 
and  plotted,  a  new  temperature  estimate  will  give  points 
B.  It  is  easy  to  see  how  they  will  converge  to  a  solution. 

The  intersection  of  the  lines  on  the  negative  slope  of  the 


Pic.  16  -  Speed  vs.  temperature  -  converging 


momentum- cor. t  mu i ty  lire  is  a  physically  impossible  solu¬ 
tion  (flow  goes  from  sub  to  supersonic)  and  is  not  looked 
into  in  this  report. 

When  ioeufc  aJUibtow  lea  Us  fo  choketf  flow  (ini 'tut 
Cj  a  »  ±  joe  ect  i  s  too  high))  there  i  >  >}  c  solution  as 

Shown  k7 /’^./7.  i”vi  this  ccise.  for  the  iowzr 

temperature  estimates  the  speeds  are  obtainable  but  never 
converge.  For  higher  temperature  estimates  negative  square 
roots  are  obtained  when  running  the  program.  Foints  A,B, 
and  C  in  fig.  17  show  the  first  three  temperature  estimates. 
After  a  new  estimate  from  points  C,  the  equation  for  the 
speed  derived  from  energy  soon  produces  this  negative  root 
due  to  the  temperature  increasing  while  the  maximum  value 
allowable  from  energy  has  already  been  reached.  A  glance 
at  the  energy  equation  in  the  next  subsection  will  illustrate 
this  further.  If,  this  occurs  then  the  speed  of  the  culv* 
or  fuel  flow  must  be  decreased  to  obtain  a  solution. 

3.  Program  Analysis 

The  following  detailed  discussion  of  the  program  will 
enable  the  reader  to  fully  comprehend  the  technique  used 
in  solving  this  problem.  Each  step  will  be  numbered  for 


reference.  The  reduced  sensible  enthalpies, 


/H-E.f 
\*T  k 


T 

formation  enthalpies,  /hf\  ,  dimensionless  specific  heats, 

\<*/i 

T  T 

jCDj  ,  and  entropy  coefficients,  [Cs|  ^  which  are  needed  can 


be  found  in  thermodynamic  tables  (  JANNAF,  AFCSR  TR  3641) 


•  a*.  • \X. 


for  temperatures  ranging  from  100  to  6000  K.  A  inter- 

A 

polation  method  is  usecl  to  calculate  values  for 
temperatures  between  those  in  the  tables.  In  general  the 
program  leads  to  a  calculated  throat  area.  This  area  is 
compared  to  a  specified  area  and  the  upstream  Mach  #  is 
modified  until  the  calculated  area  is  equal  to  the  speci¬ 
fied  area. 

The  program  begins  by  establishing  a  steady  combustion 
for  a  fixed  fuel  flow.  The  areas  of  the  primary  flow  and 
fuel  line  are  known.  The  mass  flow  of  the  air,  as  well  as 
as  the  stagnation  temperature  and  pressure  of  both  the  fuel 
and  the  air  before  they  are  mixed  are  specified.  Since  the 
injection  speed  of  the  gaseous  fuel  is  sonic,  fuel  flow 
changes  occur  from  changes  in  the  fuels'  stagnation  pres¬ 
sure.  Gaseous  hydrogen  has  been  chosen  as  the  fuel. 


M  .  is  calculated  before  the  fuel  is  added. 
c,l 


1:  M_ 


m. 


JyR alrTc , l7 ( Pc , lAc , 1Y } 


‘c,l  “‘air 
The  mass  flow  of  the  fuel  is  found. 


2:  mf  ■  AfPf/  lRfTf 


The  next  step  is  skipped  the  first  run  through. 


•  o  _  r  —  a  ~  Y+l 

3=  "air  *  Ac,lp=,l^  Mc,l/(  Mc,1)7T7-1 


)  ) 


The  fuel-air  flow  mixture  ratio  is  then 


mr  ‘“air 


Next  the  static  pressure  of  the  fuel  as  well  as  the  static 
temperature  and  pressure,  and  the  speed  of  the  upstream 
primary  flow  (air)  calculated. 

0  ,  Yf/Yf-1 

5:  Pf  -  Pf  /  (1+1-1  M*)  r;  1 


6=  Tc,l  -  Tc.l  / 


7:  Fc  1  *  K  1  (To  ,)Y/ir"1 

C,1  C  ,  1  C  ,  1  C  ,  i 


8  =  uc,l  ■  Mo,l  JV1'1 


The  maximum  speeds  possible  for  the  combustion  gas  are 
calculated.  They  are  based  on  momentum-continuity  and 
on  energy  and  lead  to  the  two  speeds  discussed  earlier  in 
subsection  C. 


9:  U 


(M-C) 
c,2  max 

(E) 


(f(Uf+RfTf)+(Uc>1+RalrTCtl))  /  (1+f) 

“*T 


- - - - - 

10:  ^  o  mav  »  (aUR^-CH-EoTf  +U*)f)+  A)/  (1+f) 
c,2  max  ^  f  f  fuel  ^ 


where 


*  *  Vm 

u  A  air  2 


Now  the  moles  of  Oj*  air  and  the  moles  of  the  combustion 
gas  are  found.  These  will  lead  to  the  value  of  the 
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>  \  • 


formation  enthalpy  of  the  combustion  gas 


U:  v0j  -1Hf  /  (4.76  f»valr) 


12:  v&lr  •  /  <f*alr> 


13 :  V  *  ''air  +  1/2 


The  gas  constant  of  the  combustion  gas  is  then 


14:  Rcg  *  HalrO-+^^alr/2TAf))/(l+n 


The  formation  enthalpy  per  unit  mass  of  the  combustion 
gas  is  given  in  terms  of  the  absolute  formation  enthalpies 
per  kmole.  These  absolute  formation  enthalpies  for  differ¬ 
ent  species  are  listed  in  table  2. 


Calculation  of  the  two  speeds  based  on  the  equations  of 
motion  is  now  possible.  The  temperature  is  iterated  on 
until  the  values  of  the  two  speeds  agree  to  a  desired 
accuracy.  A  new  estimate  of  the  combustion  gas  temper¬ 
ature  is  obtained  by  using  a  multiplying  factor  times 
the  difference  in  the  two  speeds.  This  factor  depends 
on  the  rates  at  which  the  speeds  change.  For  this  step 


a  value  of  X  *  0.03  was  used. 


The  pressure  of  the  combustion  gas  is  then  found  from 
the  continuity  equation. 


(1+f)R=gTc,2 


airTc,l  +  f  W 


C,1  C,1 


V1? 


The  conditions  at  the  point  of  choking,  the  throat, 
are  now  found.  An  estimated  temperature  for  the  throat 
is  needed  to  get  the  formation  enthalpy  per  kmole  of  the 
gas  there.  Since  the  speed  is  sonic  the  temperature  will 
be  lower  than  that  of  the  combustion  zone.  Therefore 


and  the  specific  heat  is 


2A:  “*f(1+f)  /  ( |  V 

^tair 


Since  the  stagnation  enthalpy  of  the  flow  does  not 
change,  this  value  is  calculated  at  the  throat  and  is  com¬ 
pared  to  the  value  from  the  combustion  zone.  A  new  throat 
temperature  is  estimated  based  on  the  old  value  and  the 
difference  in  the  stagnation  enthalpiy  values.  The  temper¬ 
ature  is  iterated  on  until  the  two  enthalpy  values  are 
approximately  the  same.  A  multiplying  factor  X=0.3  was 
used  for  this  iteration. 


25: 


,h  x  c,2  calc 

nn 


The 


26:  A/hf\  0,2 


VC/ eg 


,2^f\ 


Whs 


u 

^c,2  Am  +  U‘  »  /h  \  C*^  calc 

^Icg  +  V 2  (M 


T 


\5UC6 


27:  IP  ABS  (A/h 


Lc ,  2  ■ 


r)  eg 
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GO  TO  #29 


IS 


.  T4  -  I* 
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CO  TO  *21 


V •*  »  *  »  ’  »  *  *  *  s  •  •  *  .  *  i  **  »  *  »  k  «*.  V. 


The  entropy  coefficients  of  the  gas  at  the  combustion 
zone  and  at  the  throat  are  then  found  and  are  used  to 


determine  the  static  pressure  of  the  gas  at  the  throat. 

a9!  N-2  ■  ((°’)hc/ ♦  ^4)Ta;;2 - k^)oT)/Ve 

3o:  'Ml)'*** 


31:  P+ 


Pc,2Tt 


T 

(?-U 


c,2 


eg 


Now  a  calculated  throat  area  is  found.  This  is  then 

compared  to  the  specified  area.  If  the  two  are  different 

a  new  value  for  the  Mach  #  of  the  upstream  primary  flow 

(M  . )  is  calculated  and  the  entire  program  is  run  again 
c 

iterating  on  M  •  until  the  two  areas  are  the  same.  The 
c ,  x 

multiplying  factor  here  is  X®500. 


calc  #  ,  I 

32:  At  -  (mf+mair)  \J 

p 


(r««Vt,  «> 


33:  IP  ABS  (A4 


spec  calc 


-  At  )  <  10 


-10 


spec  calc 

3**  Mc,l  *  Mc,l+  X(At  -  At  > 
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GO  TO  #35 


GO  TO  #3 


•  «*  »  - 
%****••  L1 

•  %  v 


.  .‘V 


•  .  •  #  *  <*  •  % 
/  •/  /  •; 

cV-X*' 

%  *  v  • 

*•  .*•  .*»* 


i  -  •  • »  .  «  _ 


J.V.  ?T7~TT*?  '.‘  .1  ■y.iv  -«.• 


’T  *  .'  '.’  '  -•  ?  ^  r  ^  r  t  -.’  7 


k* 


The  last  calculations  give  the  stagnation  pressure 
and  temperature  of  the  combustion  gas.  An  estimate  for 
the  stagnation  temperature  is  obtained  assuming  a  perfect 
gas.  The  temperature  is  iterated  on  until  the  stagnation 
enthalpy  is  correct. 


est 

35:  Tc,2*To,2  +  !Li«i 


eg 


36:  /h 


/hf> 


‘c,2  calc 


eg 


'Hf\  °’2  *  ^alr^A  °’2  * 


WJ 


air 


ifHr 


\/v 


eg 


T® 


37:  /h-\  °*2  calc,  -  /  \  c*2  cal° 

(r)  *  Te.2(M  ^cg 

\^/cg  W/cg 


T  °  T  T° 

38:  A  /h^N  c*2  -  /h^v  c»2  -  /hf\  c»2  calc 

W  CE  MU  eg  VF7  eg 


The  first  value  on  the  right  side  of  eq.  38  is  obtained 
from  line  25. 


39:  IF  ABS  (A/hjN  c*2  ) 

\C/cg 


<  10 


-7 


GO  TO  #l»  1 


t0!  To,2  ’  Tc,2  +  °'2) 

VlC'/cg 
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Finally 


43:  PRINT  VALUES 
M4:  STOP 


Figures  18-31  show  the  results  of  a  program  run 
using  the  equations  listed  on  the  previous  pages.  The 
fuel  used  was  gaseous  hydrogen. 


4.  Results 

Figures  18-31  illustrate  the  effect  of  a  fuel  flow 
p*op**li*s  t*bi»*i* 

change  on  -teupstrea  w  ^  flow.  It  is  evident  that  the  up¬ 
stream  conditions  change  due  to  the  downstream  addition  beat  (vkoL 
vfutr  caM 

y»vv*s.Even  though  the  pressure  rise  was  small  the 

speed  decreased  by  20?  and  the  mass  flow  by  20?. 

Presently  an  experimental  device 
is  being  constructed  to  these  results. 
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'able  2  -  Thermodynamic  Values  of  Air 


Temp. 

T 

(M 

& 

\  1 

yFF/air 

yft  f  air 

100 

4.24293 

3.48805 

3.50021 

200 

4.14739 

3.49302 

3.50142 

300 

4.10166 

3.49805 

3.50951 

400 

4.07375 

3.50441 

3.53964 

500 

4.05556 

3.51745 

3.59503 

600 

4.04349 

3.53555 

3.67121 

700 

4.03564 

3.56180 

3.75530 

800 

4.03084 

3.59063 

3.83840 

900 

4.02821 

3.62299 

3.91735 

1000 

4.02701 

3.65572 

3.98804 

1100 

4.02689 

3.68882 

4.05031 

1200 

4.02758 

3.72108 

4.10616 

1300 

4.02870 

3.75313 

4.15385 

1400 

4.03017 

3.78297 

4.19756 

1500 

4.03199 

3.81218 

4.23507 

1600 

4.03384 

3.83960 

4.26844 

1700 

4.03593 

3.86581 

4.29844 

1800 

4.03801 

3.89023 

4.32586 

1900 

4.04025 

3.91302 

4.35070 

2000 

4.04236 

3.93666 

4.37358 

2000 


Table  3 

Absolute  Formation  Enthalpies 
For  Eleven  Species 
At  0  Degrees  Kelvin 
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